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^P" 1 ' Abstract 

It is argued that within the correct treatment of analytical properties of the transition 
amplitudes, in particular, the second order pole structure, characteristic for the n — n tran- 
sition in nuclei, the "infrared divergences" discussed in some papers, do not appear. Ex- 
plicit calculation with the help of diagram technique shows that the neutron-antineutron 
I/-) 1 oscillations are strongly suppressed within a deuteron, as well as within arbitrary nucleus, 

in comparison with the oscillations in vacuum. 

in 

1. The neutron-antineutron transition induced by the baryon number 
violating interaction (AB = 2) predicted within some variants of grand unified 
theories (GUT) has been discussed in many papers since 1970 [1], see [2, 3, 
4, 5, 6].. Experimental results of searches for such transition are available, in 
vacuum (reactor experiments [7], and references therein), in nucleus 16 [8] and 
in Fe nucleus [9], see also the PDG tables. 

During the later time there have been many speculations that the neutron- 
antineutron oscillations in nuclei are not suppressed in comparison with the n—n 
transition in vacuum [10, 11]. The arguments were based on the "true field- 
theoretical approach" to this problem. The result of [10] has been criticized in 
a number of papers [12, 13, 14, 15] which used somewhat different approaches 
(potential, S-matrix, diagram), and general physics arguments. 

However, in view of continuing publications [11] containing same state- 
ment as in [10], it seems to be necessary to analyze this problem just within the 
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quantum field theory based approach used in [10, 11]. My consideration is close 
to the approach of paper [14] where the diagram technique has been applied 
to study neutron- antineutron transition in nuclei, although differs from [14] 
in some details. More recent realistic calculations of the neutron-antineutron 
tansition in nuclei can be found in [16, 17]. 

Here we give first some general arguments based on analytical properties 
of amplitudes in favour of suppression of n — n transition in nuclei (section 3). 
The simplest example of the deuteron where the final result can be obtained in 
closed form, is considered in details in section 4, and the result of [14] for the 
case of the deuteron is reproduced. The analogy with the nucleus formfactor 
at zero momentum transfer is noted in section 5. 



2. To introduce notations, let us consider first the nn transition in vacuum 
which is described by the baryon number violating interaction (see, e.g. [2, 13, 
14]) V = n n n(Ji/2, o\ being the Pauli matrix. \i nfl is the parameter which has 
the dimension of mass, to be predicted by grand unified theories and to be 
defined experimentally 1 . The n — n state is described by 2-component spinor 
lower component being the starting neutron, the upper one - the appearing 
antineutron. The evolution equation is 

= (Vo + V)t (1) 

with Vo = rriN — i~f n /2 in the rest frame of the neutron (m^v is the nucleon 
mass, 7„ - the (anti)neutron normal weak interaction decay width, and we take 
In = In, as it follows from CP-invariance of weak interactions). Eq. (1) has 
solution 

V(t) = exp [-i (iJ nn tai/2 + V t)] * = 

(2) 

Here is the starting wave function, e.g. \I/o = (0, 1) T . In this case we have 
for an arbitrary time 

^(n, t) = —isin^^-exp(-iVot), ^(n, t) = cos ^""^ exp(-iVot) , (3) 

which describes oscillation n — n. Evidently, for large enough observation times, 
t obs ^> I /(i nn , the average probabilities to observe neutron and antineutron are 

1 There is relation fj, nn = 25m with the parameter 8m introduced in [2]. The neutron-antineutron oscillation 
time in vacuum is r nft = l/5m — 2//x nS , see also [17] and references in this paper. 
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equal if we neglect the natural decay of the neutron (antineutron): 



W(n) = |*(n)| 2 = |*(n)| 2 = W(n) = 1/2. (4) 

This case is, however, of academic interest, only, since j n ^> \i nn 2 It should be 
stressed that in the vacuum neutron goes over into antineutron, also discrete 
localized in space state, which can go over again to the neutron, so the oscillation 
neutron to antineutron and back takes place. 

Since the parameter \i nn is small, expansion of sin and cos can be made 
in Eq. (3) at not too large times. In this case the average (over the time 
t obs <C 1/ Unn) change of the probability of appearance of antineutron in vacuum 
is (for the sake of brevity we do not take into account the (anti)neutron natural 
instability which has obvious consequences) 

..2 fobs 

W(n-t obs )/t obs = \^(n,t obs )\ 2 /t obs ~ (5) 

which has, obviously, dimension of the width T. So, in vacuum the transition 
n — > n is suppressed if the observation time is small, t obs <C l/// n n- From 
existing data obtained with free neutrons from reactor the oscillation time is 
greater than 0.86 10 8 sec ~ 2.7 years [7], therefore, 

li nn < 1.5KT 23 eV, (6) 

very small quantity. 

Recalculation of the quantity \i nn or r nn from existing data on nuclei 
stability [8, 9] is somewhat model dependent, and different authors obtained 
somewhat different results, within about 1 order of magnitude, see e.g. discus- 
sion in [14, 16, 17]. Most recent results for fi nn obtained from nuclear stability 
data are close to (6) [16, 17], see also next section. 

3. In the case of nuclei the n — n line with the transition amplitude 
linn is the element of any amplitude describing the nucleus decay A — > (A — 
2) + mesons, where (A — 2) denotes a nucleus or some system of baryons 

2 It is a matter of simple algebra to calculate the integrals over time of the probabilities \^>(n, i)| 2 and 
l*(R,()l 2 :. 

f w "' i,|2dt = 5^OT f '•<»•«>''* "5itf&EJ 

for neutron as initial state and for arbitrary, but different from zero j n . The difference between both quantities 
is obvious, and disappears when j n — > 0. 
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with baryonic number A — 2, see Fig. 1. The decay probability is therefore 
proportional to /j^ n , and we can write by dimension arguments 



V(A ->> (A - 2) + mesons) 



m ' 



(7) 



where mo is some energy (mass) scale. For the result of [10, 11] to be correct 
the mass mo should be very small, mo ~ fi n n ~ 10~ 23 eV, but we shall argue 
that mo is of the order of normal hadronic or nuclear scale, mo ~ rrihadr ~ 
(10 — 100) MeV. We can obtain the same result from the above vacuum formula 
(5), if we take the time t obs ~ 1 /rrihadr ■ 




mesons 



Figure 1: The diagram describing n — n oscillation in nucleus A with subsequent annihilation 
of antineutron to mesons. The final state has the baryon number A — 2. 



Indeed, the matrix element of any diagram containing such transition 

T(A ->> (A - 2) + mesons) ~ 

^ r T ,, , / - Tin + (A — 1) — > (A — 2) + mesons) ,^ 
~ - Z) / V(^; n, (A - ^ _ £0 + J)2 ~ dE n - 

^- 2m (A-Z) d -^(E n = El), (8) 

according to the Cauchy theorem known from the theory of functions of complex 
variable. E n is the neutron (antineutron) energy - integration variable, E® is the 
( anti) neutron on-mass-shell energy E® ~ m^^y 1 j2m^. The energy-momentum 
conservation should be taken into account for the vertex V(A — > n + (A — 1)) 
which includes the propagator of the (A — 1) system, and for the annihilation 
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amplitude T. The case of the deuteron considered below is quite transparent 
and illustrative. 

The amplitude T which describes the annihilation of the antineutron, and 
the vertex function V are of normal hadronic or nuclear scale and cannot, in 
principle, contain a very small factors in denominator (or very large factors, of 
the order of 10 15 , in the numerator). By this reason we come to the above Eq. 
(7), and the resulting decay width of the nucleus is very small, 

T{A ->• (A — 2) + mesons) < 10~ 3 Vnn, (9) 

at least 30 orders of magnitude smaller than the inverse time of neutron- 
antineutron oscillation in vacuum \i nfl . From Eq. (7) or (9) we obtain 

fi nfi ~ \JT\A — > (A — 2) + mesons)mo, (10) 

and when one tries to get restriction on \i nfl from data on nuclei stability [8, 
9] the result is close to that from vacuum experiment [7], somewhat smaller, 
within one order of magnitude [5, 6, 14]. The result of [16] differs from that 
of [14] for heavier nuclei, and the authors [16] come to the conclusion, that 
experiments with free neutrons from reactor could provide stronger restriction 
on the neutron-antineutron transition parameter than experiments on stability 
of nuclear matter . 

According to [10, 11] the probability of the nucleus decay is proportional 
to W(t obs ) (t obs ) (the process proceeds similar to the vacuum case), 

where t obs is the large observation time, of the order of ~ 1 year or greater. By 
this reason the extracted value of \i nfl is smaller than that given by Eq. (10), 
by about 15 orders of magnitude. Technical reason for strange result obtained 
in [10, 11] is the wrong interpretation of the second order pole structure of 
any amplitude containing the n — n transition. Instead of using the well de- 
veloped Feynman diagram technique, the author [10, 11] tries to construct the 
space-time picture of the process by analogy with the vacuum case, which is 
misleading, see also discussion in conclusions. 



3 Thcre is, in fact some kind of competition between both methods, and final result will depend on the 
progress to be reached in both branches of experiments — with free neutrons and with neutrons bound in 
nuclei. Friedman and Gal [17] obtained the restriction r„ ft > 3.310 8 sec from the latest datum on le O stability 
and using the potential approach. Experiments with ultracold neutrons in a trap have been proposed and 
discussed in [3, 18], but not performed till now. 
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4. We continue our consideration with the case of the deuteron which 
is quite simple and instructive, and can be treated using standard diagram 
technique 4 . The point is that in this case there is no final state containing 
antineutron - - it could be only the pn state, by charge conservation. But this 
state is forbidden by energy conservation, since the deuteron mass is smaller 
than the sum of masses of the proton and antineutron. Therefore, if the n — n 
transition took place within the deuteron, the final state could be only some 
amount of mesons. The amplitude of the process is described by the diagrams 



d 




mesons 



Figure 2: The diagram describing n — n oscillation in the deuteron with subsequent annihilation 
of antineutron and proton to mesons. 



of the type shown in Fig. 2 and is equal to 

, x r T(np — >• mesons) d A p . s 

T(d -f mesons) = lgi .„ m j {f _ mj , )[(d _ p)2 _ m%]2 J^y- (". 

The constant gd np is normalized by the condition [20, 21, 22] 



2 

9dnp ^ 



lo7r mjsi \ m^ 

which follows, e.g. from the deuteron charge formfactor normalization Fd(t = 
0) = 1, see next section, k, = ^/m/ved, Cd — 2.22 MeV being the binding energy 
of the deuteron. For the vertex d — > np we are writing Zm^gdnp to ensure the 
correct dimension of the whole amplitude. 

The integration over internal 4- momentum d A p in (11) can be made easily 
taking into account the nearest singularities in the energy p$ = E, in the non- 
relativistic approximation for nucleons. As we shall see right now, the integral 

4 It has been considered in fact in [19] within reasonable framework of diagram technique, however, the author 
has drawn later wrong conclusions from this consideration. 
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over d?p converges at small p ~ k which corresponds to large distances, r ~ 1/k. 
By this reason the annihilation amplitude can be taken out of the integration 
in some average point, and we obtain the approximate equality 

T(d — > mesons) = gdnp^N^nn^dNNT^np — > mesons) (11a) 

with 



0? 4 J9 

' ,/V V : J^y-I (p 2 - m 2 N )[(d - p) 2 - m%} 2 



d A p 



(27r) 4 (2m) 3 ■/ (p — m N — p 2 / (2m N ) + iS)(md — m N — p — p 2 /(2m N ) — iS) 2 

= I (1S) 

(2irf8m N [K 2 + p 2 } 2 64™^' 

This integral converges at small \p\ ~ k, more details can be found in the next 
section. The decay width (probability) is, by standard technique, 

T(d — > mesons) ~ lAin9dnp^dNN m N j \T(np — > mesons)\ 2 d§ (mesons), (14) 

^(mesons) is the final states phase space. Our final result for the width of the 
deuteron decay into mesons is 

^mesons ^ ^< N^WL^O ^ 3Z7>* \Pcm.°%\ m .-+0 , (15) 



where p c . m . is the (anti)nucleon momentum in the center of mass system. This 
result is very close to that obtained by L.Kondratyuk (Eq. (17) in [14]) in 
somewhat different way, using the induced np wave function 5 . 

The annihilation cross section of antineutron with velocity vq on the pro- 
ton at rest equals 

1 r 

a(np — > mesons) = 2 — / \T(np — V mesons)] 2 d& (mesons). (16) 
According to PDG at small vq, roughly, Vqct^ 71 v — (50 — 55)m6 ~ (130 — 



140) GeV 2 . So, we obtain fi nn < 2.5 10 24 eV, or r nfl > 510 8 sec if we take 



3 The result Eq. (17) in [14] can be rewritten in our notations as 



r^„ s - 0.01 M L— [vo<r ann np] v . , (17') 
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K 



which differs from our result by some numerical factor, close to 1 and not essential for our conclusions. 



optimistically same restriction for the deuteron stability as it was obtained 
for the Fe nucleus, Td ~ T Fe > 6.510 31 yr [9]. Our result (15) is valid up 
to numerical factor of the order ~ 1, since we did not consider explicitly the 
spin dependence of the annihilation cross section and the spin structure of the 
incident nucleus. Same holds in fact for the results obtained in preceeding 
papers, see [2, 14] e.g. 

Additional suppression factor in comparison with the case of a free neu- 
tron is of the order of 

jlnn/Ki ~ 1CT 31 

in agreement with our former rough estimate (9), and disappears, indeed, when 
the binding energy becomes zero 6 . The binding energy of the deuteron should 
be very small, to give n ~ fj, n n, to avoid such suppression. At such vanishing 
binding energy the nucleons inside the deuteron are mostly outside of the range 
of nuclear forces, similar to the vacuum case. 

Results similar to (15) can be obtained for heavier nuclei, see [5, 6, 14, 16, 
17]. The physical reason of such suppression is quite transparent and has been 
discussed in the literature long ago (see [2, 13, 15] e.g.): it is the localization 
of the neutron inside the nucleus, whereas no localization takes place in the 
vacuum case. In the case of the deuteron or heavier nucleus the annihilation 
of antineutron takes place, and final state is some continuum state containing 
mesons. By this reason the transition of the final state back to the incident 
nucleus is not possible in principle, and there cannot be oscillation of the type, 
e.g. d — > mesons — > d. This is important difference from the case of the free 
neutron. 



5. As we noted previously, the presence of the second order pole in 
intermediate energy variable is characteristic for the processes with the neutron 
- antineutron transition, but it is in fact not a new peculiarity, it takes place also 
for the case of the nucleus formfactor with zero momentum transfer, Fa{q = 0). 
Let us consider as an example the deuteron charge formfactor. In the zero range 
approximation it can be written as 



i(2mg dnp ) 2 r d 4 p 



Fd ^ (2tt) 4 J (p 2 - ml) [(d - pf - m%] [(d-p + qf - m%] ' 



6 There is no final formula, for Td^mesons in [19] to be compared with our result (14), (15). Numerically, 
however, the result of [19] is in rough agreement with our and [14] estimate. 
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Behind the zero range approximation g dnp should be considered as a function of 
the relative n — p momentum, not as a constant. For q = second order pole 
appears, and we come to the expression for F(q = 0) containing the integral 
IdNN introduced above in Eq. (13): 

F d (0) = (2m N g dnp ) 2 I dNN . (18) 




Figure 3: The diagram describing the deuteron charge formfactor. 



In the nonrelativistic approximation, when only the nearest in energy 
E = po singularities are taken into account, the integral over the energy has the 
structure 

f —2iri , 

tv* ~ J (E-a + isw-b-isy = J^W (19) 

a = rriN + p 1 /2rriN. b = rrid — ttin — p 2 /2rriN, a — b = e d + p 1 /mjsr, and can 
be calculated using the lower contour, or the upper contour, with the help of 
formulas known from the theory of functions of complex variables. After this 
we obtain 

Since F d (0) = 1, this leads to the above mentioned normalization condition 

This relation between the constant g dnp and the binding energy of the 
weakly bound system (deuteron in our case) is known for a long time [20, 21, 22]. 



7 As it is known from nonrelativistic diagram technique, the wave function of the deuteron in momentum 
representation is ^d(p) = gdnp/i^ 3 ^ 2 ^ 2 +F 2 )]: therefore, the normalization of the charge formfactor F^(0) = 1 
follows from the normalization of the deuteron wave function, which is also well known from quantum mechanics. 
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It was obtained in [20, 21, 22] using different methods, dispersion relation, for 
example. 

If the infrared divergence discussed in [10, 11] takes place for the process 
of n — n transition in nucleus, it should take place also for the nucleus form- 
factor at zero momentum transfer. But it is well known not to be the case. 
There is no "new limit on neutron - antineutron transition" [10]; instead, one 
should treat correctly singularities of the transition amplitudes in the complex 
energy plane, in particular, the second order pole contribution to the transition 
amplitudes. The author of [10, 11] tries to reconstruct the space-time picture 
of the process, but the correspondence of this picture to the well justified am- 
plitude, as it appears from Feynman diagrams, is questionable. The infrared 
divergence discussed in [10, 11] is an artefact of this inadequate space-time 
picture of the whole process of n — n transition with subsequent antineutron 
annihilation. Another quite unrealistic consequence of this space-time picture 
is the nonexponential law of the nucleus decay. 

Field-theoretical description of nuclear reactions and processes is poten- 
tially useful, it allows to study some effects which is not possible in principle 
to study in other way, e.g. relativistic corrections to different observables. One 
should be, however, very careful to treat adequately analytical properties of 
contributing amplitudes. E.g., in the case of the parity violating amplitude 
of np — > d'y capture it was necessary to take into account contributions of all 
singularities (poles) of the amplitude in the complex energy plane, not only 
contributions of the nearest poles in energy variable, as it is made usually in 
nonrelativistic calculations. The nonrelativistic diagram technique developed 
up to that time turned out to be misleading for the case of physics problem 
considered in [23]. Besides, and it is the spesifics of the processes with photon 
emission, the contact terms should be reconstructed to ensure the gauge invari- 
ance of the whole amplitude of photon radiation [23] . 
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